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Abstract 

Recently a new class of time-dependent Bell inequalities in Wigner form was introduced. The 
structure of the inequalities allows experimental studies of quantum and open quantum systems 
in external fields. In this paper we study the properties of the time-dependent Wigner inequalities 
using the time evolution of neutral pseudoscalar mesons. It is shown that it is always possible to 
find a range of parameters to test for violation in an experimentally accessible area. The effect of 
the relaxation of the inequalities for large time scales is demonstrated. 
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I. INTRODUCTION 


In 


0 


a new class of time-dependent Bell inequalities in Wigner form was introduced. The 
structure of the inequalities allows experimental studies of the quantum and open quantum 
systems in an external fields. In this paper we study the properties of the time-dependent 
Wigner inequalities using the time evolution of neutral pseudoscalar mesons. 

In Q the following question was raised for the first time: can the properties of a macro¬ 
system which in quantum theory are described by noncommuting operators be simultane¬ 
ously the elements of reality (i.e. to exist simultaneously), even if these properties can not be 
measured by any macro-device? According to the Copenhagen interpretation of quantum 
mechanics, the answer is no Q|. Bohr demonstrated the fundamental difference between 
the statistical and the Copenhagen interpetations, but did not give a conclusive proof [2:]. 


An attempt to move the problem of the simultaneous existence of the elements^) 
reality from the gedanken to the experimental realm has been made by J. Bell 


physical 
6|. Bell’s 


idea has been further developed by Clauser, Horne, Shimony, and Holt 7]. Since then the 
idea has been thoroughly studied, developed, and criticized by many Bell’s proponents and 
opponents B-Q. 


We address here the question of how to express the fact, that some set of a micro-system 
characteristics (observables) is simultaneously the set of elements of the physical reality, 
even if that set can not be simultaneously measured by any macro-device. One possibility is 
to assert that the joint probability of simultaneous existence of members of the set is non¬ 
negative. In quantum mechanics the probability is the universal object. Such a proposition 
(however without the clear statement) was used by Bell in { 4 ]. The idea of the non-negativity 


of the joint probability was proposed by E. Wigner 


12J- 


Following Wigner’s approach, let us suppose that a quantum system decays at the time 
t 0 into two subsystems “1” and “2”, each having three observables a, 6, and c. Let each 
observable be dichotomic variable (able to have only two distinct values). For simplification 
let us set these values to ±1. We will use the following definitions: if the observable a is 
equal to +1 we denote this as a+\ and so on. At the time t 0 let all three observables to 
satisfy the anticorrelation condition 

n±\t 0 ) = — (t 0 ), (1) 

where = {aW, b^\ c^} and i = {1, 2}. 
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An example of such a system is a pseudoscalar particle which decays into a fermion- 
antifermion pair with a Hamiltonian 

n (PS \x) = g <p(x) (f(x)j 5 f(x)) N , (2) 

which automatically provides full anticorrelation of the fermions’ spin projections onto any 
direction. In fl2]), the <p(x) is a pseudoscalar held, and f(x) and f{x) are fermionic helds. 

Let a, b , and c exist simultaneously as the elements of physical reality, i.e. any of their 
double and triple joint probabilities are non-negative. Then, under the assumption of locality 
at time t 0 and using Kolmogorov’s axiomatics of probability theory, the following inequality 
can be obtained 111: 


w o) < w ^cV 2) ,&+ 1} ,to) + w (a^\c^\t 0 ^ . 


If one drops the t 0 from (J3J), it transforms into the Wigner inequality [12] 

w (a+\b+^ < w (c+\b+^ + w . 


(3) 

(4) 


We name the inequalities ([3]) and ]4|) time-independent or static Wigner inequalities, to 
distinguish them from the time-dependent inequalities in [] . The observables a, b, and c may 
correspond to non-commuting operators and, hence, cannot be simultaneously measured by 
any macro-device. As was mentioned in 13], the Wigner inequalities are more suitable for 
experimental tests, due to the fact that the probabilities, unlike the correlators, are well 
defined in both non-relativistic quantum theory and in quantum held theory. 


In 


0 . 


a new class of Wigner inequalities was obtained with a direct dependence on time: 


w (a { l\b { l\t^ < ( 5 ) 

< w (a+\to) ->• (w (b^\to) ->• b+\t)j + w ^ U) (t 0 ) ->• b^\t)^ w c+\ t 0 ) + 

+ w (aL 2) (t 0 ) ->• a+ } (^)) (b+\t 0 ) ->• b { +\t)^ + w (V_ 1} (t 0 ) ->• w (a-\c+\t > 0 ) + 

+ w (b^\t 0 ) ->■ b ( +\t)^ (w (a^\t 0 ) ->• + w (aL 2) (t 0 ) ->• a+ 2) (^))) w (c^\b^:\t 0 ) + 

+ w (b [ l\t 0 ) ->■ b+\t)J (w (a+\t 0 ) ->• a+ } (t)) + w (a^\t 0 ) ->■ w (c+ } , 6 (i) , t 0 ) • 

For closed quantum systems, w ^a^(f 0 ) —>■ = w (b^\t 0 ) ► b+\t)^ = 0, while 

w (ci+\t 0 ) — > a+\t)j = w (b+\t 0 ) —> b+\t)^ = 1. Hence ([5]) reduces to Q, as it should be 
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from physical point of view. The inequality (l3|) . in turn, is equivalent to the time-independent 
inequality ([4]). 

The Leggett-Garg —ties are based on the idea of —opic teahsat Q. They 
resemble Bell inequalities [7j but, instead of the simultaneous correlation of the two observ¬ 
ables, they concern the correlation between the values of a single observable at different 
points of time. In experiment they are closely related to the weak (non-invasive) measure¬ 
ments 115, 161, involving, for example, nanomechanical resonators 1171 or discrete lattices 

H. 

The main distinction between the Leggett-Garg inequalities and (J5]) is the fact that 
the test of (J5J) does not require weak measurements. The measurement is fully invasive. 
That opens the experimental possibility to verify (J5J) at contemporary high energy physics 
detectors like LHCb, ATLAS, CMS, and Belle II. 

We study here the violation of the time-dependent inequality (|5|) in systems of neutral 


pseudoscalar mesons. A 


Wigner inequalities 


fl 7 


any papers attempt to include time-dependence into the static 


12], and studies of the obtained time-dependent 


quantum theory are available 19] - 30]. A number of authors 22] - 26] try to adapt 


inequalities in 


Wigner inequalities for oscillations of pseudoscalar mesons, usually the neutral /L-mesons. 
First, these adapted inequalities are studied in terms of “flavour”-“CP-violation”-“states 


with defined masses and lifetimes”. This idea was introduced in 


22], The time-dependence is 


included by substitution of probabilities calculated in the framework of quantum mechanics. 
In this case Wigner inequalities become inequalities among the parameters e and e' of CP - 
violation. The violation of these inequalities is small and is currently beyond experimental 
reach 22], Secondly, there are attempts to include additional correlation functions which 


depend on time difference 


26]. A third way, introduced in |l9j] - 2l|, is based on the 


requirements of the causality principle and locality; however the obtained inequalities are 
not general and are suitable only for the specific situation of the oscillations of neutral 


0 


mesons. Finally, in 26] - 30], special versions of time-dependent inequalities in the form 
are introduced, but there are certain difficulties with their violation in quantum mechanics. 

To demonstrate the distinction between 0 and (j4j), we will apply them to the problem 
of oscillations of neutral pseudoscalar mesons M = {A', D , B q }, q = { d , s}. In this case the 
static inequalities ® are either not violated at all, or the scale of the violation is beyond 
experimental reach 22], The violation of the time-dependent inequality (EJ), on the other 
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hand, can be significantly enchanced by a proper choice of parameters and hence allows 
experimental tests. 


II. STATIC WIGNER INEQUALITIES FOR OSCILLATIONS OF NEUTRAL 
PSEUDOSCALAR MESONS 


The key idea of the static Bell inequalities for the task at hand was suggested in 
and developed in 


25 


26 


19 


22 


231 


The essence of the idea is that there are three naturally provided “directions” whose 
projection operators do not commute. The first is the flavour of a pseudoscalar meson. For 
example, for the B q mesons, we consider the projections onto the states | B q ) = \bq) and 
| B q ) = | bq). We define the operators for charge ( C ) and spatial (P) conjugation onto the 
states in the flavour space as 


CP | M) = e ia \ M) and CP\M) = e~ ia \M), 


where the a is a non-physical arbitrary real phase of CP-violation. This phase should be 
excluded from any experimentally testable inequalities. 

The second “direction” is the states with defined values of CP- parity, i.e. the states 

|M 1 ) = -h (|M)+e”|M», M 2 ) = -h (| M) — e ia \ M )), 

which have positive and negative CP- parity accordingly. 

The third “direction” is defined by the states with fixed values of mass and lifetime 

\M L )=p(\M)+e ia -\M) ] and I M H ) = p (\ M) — e ia - \ M) 

V v ) \ v 

The latter two states are the proper vectors of the non-hermitian hamiltonian (for which 
CPT-symmetry is preserved): 



H 



(m — i / 2 V {mi 2 — */ 2 T 12 ) e lOL 

l (ml 2 — */2 Tj; 2 ) e ia m — */2 T 


with the proper values 


E l = m L - */2 T L = TL - a/ H 12 H 2 i =H + i/pH 12 and 
Eh = mji — */2 T h = P + \/ Bti 2 H 2 i = P — i/pH\ 2 
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accordingly (here and subsequently we use the natural system of units in which h = c = 1). 
The states | Ml ) and | Mh ) are not orthogonal to each other. The complex coefficients p 
and q are subjected to the standard normalization condition: 

(M l \M l ) = (M h \M h ) = \p\ 2 + \q\ 2 = 1. (6) 


We define 


AM = M H - M L = - 2 Re ( - H 12 

\P 

Ar = Ttf-r L = 4Im \- h 12 

\p 

Note that the definition of AT in the current work is oppositely signed relative to the 
definition in 31]. 

To automatically satisfy we introduce a new variable (3, for which 

\p\ = cos/3; |g|=sin/3; and - = tg/3 = re 1 **, (3 G [0, 7 t/2 ] . 


Then 


| M l ) — p (| M) + e ^ a M) j\p ^ and | ][f R ) = p (|M) - e 4(Q+ ^ tg/5| M )) . 

Decay of a neutral vector state 1 into an MM -pair (e.g. 0(1020) —> KK or T(4S') —y 
BB) defines a flavour-entangled wave function of the MM- system at t — t 0 : 

| T(t 0 )) = (|M) (2) |M) (1) - |M) (2) |M) (1) 


v/2 V 


(7) 


The distinction between the hrst and the second meson can be provided by their direction 


in the experimental device, see for example 


22 


25 


26]. 


We obtain the static Wigner inequalities following the logic of 26]. We make the corre¬ 
spondence a+ — > Mi, a_ — > M 2 , b + — > M, 6_ —* M, c+ —> Mh and c_ —> Ml- Then (]4]) 

becomes: 


w(M?\ M (1) , t 0 ) < w(M[ 2 \ M%\ t 0 ) + w(M {2 \ M«, t 0 ). 
Substitution of probabilities (lAljl from Appendix |A] into (JHJ) leads to: 


( 8 ) 


\q \ 2 - H 2 < \p + q I 2 - 


(9) 
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As there is no unambiguous correspondence between the projections of the meson states 
onto various “directions” and the projections from (J3J), we set b + —» M and —y M while 
keeping the a± and c±. Then (J3J becomes: 


M' 1 ', t 0 ) < w(m[ 2 \ Mj}\ to) + w(M%\ tf», t„), 


( 10 ) 


and, taking into account (lAlj) : 


\p I - M < \p + q\ 


( 11 ) 


One can merge (j9]) and (fill) as: 


\p\ - M < \p + q r 


( 12 ) 


Now let a + —> M 2 and a_ —>■ Mi, keeping the b± and c±. Then from ([4j) follows 

w(M?\ M' 1 ), to) < w{M®, M£\ t 0 ) + w{MP, M (1) , t 0 ), (13) 

and, taking into account (1A1|) : 


|g| - \p\ < | p-q\‘ 


(14) 


Finally let a + —* M 2 , a_ —> Mi, b + — > M, b_ — > M, and c + — > M H , c_ —* M L . Then, from 
(HI) it follows that 

W (Mf, M«, t 0 ) < w(M?\ m£\ to) + , t 0 ), (15) 


and, in turn, 

\p \ 2 ~ \q \ 2 <\p~ q \ 2 ■ (16) 

The inqualities (fl4l) and (fTfil) can be merged into 

\\p \ 2 - \q\ 2 \ < \p - q \ 2 ■ (17) 

The inequalities (1T21) - (1X71) represent the full set of the static Wigner inequalities for the 
oscillations of neutral mesons. The set is obtained from all possible correspondences between 
the a±, b±, and c± and the projection of the meson states onto the “directions” of flavour, 
CP, and the states with fixed masses and lifetimes. It is obvious that era and to can 
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TABLE I: 


mesons 


work and 


Experimental values of the oscillations and the CP-violation for the neutral pseudoscalar 


3l(]_. The minus sign of AT is due to the difference of the definitions between the current 


31]. The dimensionless variable A = 


Meson 

AT (MeV) 

AM (MeV) 

tan/3 = \i/ P \Z P 

A 

B° 

K° 

D° 

-6.0 x 10~ n 

-7.3 x 10~ 12 

-2.1 x 10~ n 

1.2 x 10” 8 

3.5 x 10~ 12 

-6.3 x 10 -12 

1.0039 ±0.0021 

0.99668 ± 0.00004 

0 92 +0 ' 12 

u ' y ^-0.09 

-0.2 x 10 3 

-4.8 x 10" 1 

0.3 


not be violated simultaneously - one can sum (1T2|) and (TT71) . However there are no physical 
arguments to prefer (fT2l) over (fT7l) or vice versa. 

Using (1A1I) it can be shown that (fT2l) and (1771) reduce to 

|cos2/3| ± cosC sin(2/3) < 1. (18) 

Inequality (fT8l) does not contain the unphysical phase a, as must be the case for any exper¬ 
imentally testable relation between the observables in quantum theory. 

We now check whether (fT2l) is violated in systems of neutral pseudoscalar mesons. In 
Table I the current experimental values of \ q jv\M are shown. For all the neutral mesons 
\ q /p\ e ^ p ~ 1, i.e. (3 « (3 0 = Y 4 - For the D-mesons, \ q /p\p P is not well measured, however 
within the experimental uncertainties it is consistent with one. 

For (3 = f3 0 = V 4 (be. without oscillation-induced CP-violation), (fT8j) reduces to the 
trivial inequality 

|cos Cl < 1, (19) 


which is not violated for any value of the phase C- 

However (fill does not mean that the static inequalities (ITT)) and (HT)) are never violated. 
Due to the CP-violation, the angle (3 is slightly different from (3q = Y 4 - 

For K -mesons violation of (fj~2l) was demonstrated in 22], In the case at hand, the 
coefficients p and q are defined through the CP-violation parameter e as: 


P = 


1 


1 + £ 


y/l “b l £ l 


and q = 


\/2 yU T |Y 


Then (1T2|) becomes 


Re(e:)| < 1. 


(20) 


























and (fTTj) becomes 


|Re(e)| < kl 2 , 


( 21 ) 


which corresponds to (16) from 22] if one neglects the corrections ~ |e| 2 and the modulus 


on the left hand side. The inequality 


is never violated, as |e| 


10 


-3 


31], leading 


to the upper limit |Re(e) | ~ 1CU 3 . Inequality (l2Tj) should be strongly violated, as <p e = 
(43.52 ± 0.05)° 31]. However due to the smallness of UP^violation in neutral K -mesons, 
direct experimental tests of (fl3l) and (1T51) are not possible 22]. 


III. TIME-DEPENDENT WIGNER INEQUALITIES FOR THE OSCILLATIONS 
OF NEUTRAL PSEUDOSCALAR MESONS 

We now consider time-dependent Wigner inequalities (J5]) for neutral pseudoscalar meson 
systems. Note that the normalization of probability to unity was not used in the derivation 
of ([5]). Hence ([5]) is valid for a unstable particles, whose state vector normalization is time- 
dependent. 

The time evolution of the states | Ml) \ Mh ) is given by: 

| M L (t)) = e~ iEL At | M L ) = e~ imL At ~ VL At/2 1 M L ), (22) 

| M H (t)) = e~ iEii At | M h ) = e~ imH At ~ VH At/2 1 M H ), 
where At = t — t 0 - This leads to the evolution of the states | M(t)) \ M(t) ) as: 

t |M(i)}=J+(Af)|M) - e ta ig.(At)\M) 

| ! M(t)) — 9+ (At)| M) — e~ ia -g-(At)\ M) ' 

\ q 

We then obtain the evolution of the states | M\ (t )) and | M 2 (t )) as: 

|M,(«)> = A (Y 9+ (Ai) - W(Ai)W) + e“ (g + (At) - 2j.(At)Wd . 

\M 2 {t)) = A (Y 9+ (A t) + |g_(Ai)4 | M) - e” (g+(At) + U_(Af)Wd , 

where g±(r) = - ( e~ lEHT ± e ~ lELT ). For the functions g±(r), the following is satisfied: 

\g±{r )\ 2 = e ~Y~ ( ch ± cos ( AMr )) ’ 

g* + (r)g-(r) = - ^sh A i sin (AM r)^ , 
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where F = (r# + Tl)/2. Taking into account initial condition (ITT) , it is possible to write the 
wave function of the MM -pair at arbitrary time t : 


T(t)) = e 


— p-i{. m H+m L ) At -r At 


*(*o))- 


(23) 


To simplify the subsequent calculations we from now on set to — 0, so At = t. 

To demonstrate the advantage of (J5]) over (J4]), we first consider the case of (3 = /3 0 = ^4, 
when there is no oscillation-induced CP-violation. The static inequality (0J) becomes the 
never-violated inequality (THU) . One can obtain a significant simplification by considering the 
time-dependent 
then for K— mesons, 


with the additional condition cos(j = ±1. If one neglects CP-violation, 
q\ = 1 - C 

p)k 1 + 6 


1, so cos C k = 1. For B q -mesons the effective 


Hamiltonian of the oscillations is proportional to ( V t bV. t * q ) 32], Then: 


q 

PSB a 


H- 


21 




21 


y; b v„ 

\v*v t 


tq\ 


= -l, 


hence cosCb = —1. For P-mesons the experimental results of BaBar 


and Belle 


341 


are in accordance with the assumption cos Co = 1, in which case the condition cos£ = ±1 
is well justified. 

Table II shows time-dependent Wigner inequalities for all possible correspondences be¬ 
tween the dichotomic variables a±, b ±, c± and the projections of meson states onto the 
“directions” of flavour, CP, and states with fixed masses and lifetimes. 

All calculations are performed using formulas (lAlj) and (IA2|) with the approximation 
/3 — /3o — Y 4 an d cosC = ±1. The experimental values of the oscillation parameters shown 
in Table I, and numerical estimates of the cos (m suggest the optimal choice of sets N5 and 
N6 from Table II for studying the violation of ((5]) in K- and P-mesons. For studying the 
violation of (J5]) in oscillations of neutral Bd tS -mesons one should choose sets N7 and N8. 


IV. CP-VIOLATION EFFECTS INFLUENCING THE VIOLATION OF THE 
TIME-DEPENDENT WIGNER INEQUALITIES 

We now take into account all the CP-violation effects, i.e. the case when j3 ^ /3 q, and 
cos (m ±1- Then for the various sets from Table IXTl the substitution of (IA1I) and (IA2j) 
into the time-dependent Wigner inequalities (J5j) results in eight inequalities. They can be 
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TABLE II: Time-dependent Wigner inequalities ([5]) for neutral pseudoscalar mesons with the 
approximation /3 = /3 o = V 4 an d cos£ = ±1. All possible correspondences between the dichotomic 
variables a±, b±, c± and the projections of the meson states onto the “directions” of flavour, CP, 
and the states with fixed masses and lifetimes are shown. 


N 

Correspondence of 

Time-dependent 

Violation 


the variables 

Wigner inequalities 

conditions 

1 

tt_i_ — y Mi , &_|_ — y Ad , c_|_ —y Adfj , 

1 < e" Art 

if AT > 0 


a_ — > M 2 , b_ M, c- Ml 

when cosC = — 1 


2 

—y Adi, 6 _j_ —y M , c_|_ — y Adjj , 

1 < e” Art 

if AT > 0 


a_ — » M 2 , b_ — »• M, c- Ml 

when cosC = — 1 


3 

ci-\- —^ A /25 ^ A/” , c_|_ —^ Adfj , 

1 < e” Art 

if AT > 0 


a_ —»• Mi, b- —»■ M, c- Ml 

when cost) = +1 


4 

ci-\- —^ A /25 ^ A/” 5 c_|_ —^ Mh 5 

1 < e" Art 

if AT > 0 


ci— —^ Afi, 6 — —^ A/", C— —^ Ml 

when cos£ = +1 


5 

ci-|_ —^ Mi, 6 _|_ —^ Af, c_|_ —^ Ml-, 

1 < e Art 

if AT < 0 


a_ — »■ M 2 , 6 _ — >■ M, c_ — »■ Mh 

when cos£ = +1 


6 

ci-|_ — ^ Mi, 6 _|_ — ^ Af, c_|_ — ^ Ml, 

1 < e Art 

if AT < 0 


a_ —» M 2 , » M, c_ — » Mh 

when cost) = +1 


7 

a + — > M 2 , 6 + — ■> M, c + —> Ml, 

1 < e Art 

if AT < 0 


a_ — » Mi, » M, c_ — » Mh 

when cos£ = — 1 


8 

cl- f —^ M 2 , 64 - — > AL, c + —> Ml, 

1 < e Art 

if AT < 0 


a_ — »■ Mi, b- — >• M, c- Mh 

when cos£ = — 1 



reduced to: 

1 < R n (x, r, C, A). (24) 
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The functions Rn depend on the dimensionless variables x = A Ft, A = am /af, the absolute 
value r, and the phase ( of the ratio i/p. The experimental values of AT are less or equal 
to 0 for K . 1) . and B s mesons, thus we consider only the functions R 5 - R 8 . Analytical 
expressions for these functions are given in Appendix [Bl Numerical values of the parameters 
used for the analysis of the inequalities (12T|) are given in Table I. 

We start with the system of neutral kaons. For /7-mesons the absolute values and phase 
of the CP violation parameter e are known with quite high precision. Hence r and ( are also 
well defined. In FIG. Q] the functions Rs(x, r, £, A) and Re(x, r, £, A) are shown. For kaons 
these functions are almost identical for the experimentally allowed values of r and ( (e.g. 
with r = 0.997 and ( = —0.18°, which are used in FIG |T]) . The top scale corresponds to 
the variable ct (the decay length) in mm. The bottom scale corresponds to time measured 
in units of the average kaon lifetimes z — — (T# + r^) t = Ft. Time t is calculated in the 

/7-meson rest frame. In FIG. [2] the functions R 5 . e(x, r, (, A) are shown for z < 3, i.e. in 

the most experimentally accessible area. 

Time-dependent inequalities (1241) are violated when R 5 j 6 (x, r, (, A) < 1. FIG. HJ shows 
that for the /7-mesons this violation occurs when z < 5.5, i.e., in the experimentally acces¬ 
sible area. For z > 5.5 the inequalities (1241) are not violated. This interesting effect may be 
understood if one makes an expansion of the functions Rn(x, r, (, A) by small parameters 
Ar = r — 1 and (. As an example we obtain the expansion of the function R 5 (a;', r, (, A) to 
second order. We consider Ar and ( to be of the same order of magnitude. Then: 

R 5 (x, r, C, A) » ^ {e x + 1 ) (1 - Ar) + 

+ ( 3dl2 (f) + \ Sh (I) Ch (f) “ 2 ch (|) cos ( Aa: )) ( Ar ) 2 + 

+ Qeh 2 (^j - ch (^j cos(Ax)^ ( 2 - ch Q j sin (Ax) ( Ar. 

At zeroth order in Ar and (, which corresponds to the absence of CP-violation, the effect 
of restoration of the inequality (1241) at large values of t (or z) does not appear. That is, this 
effect is fully determined by the CP-violation. A first order expansion is also not enough. 
The effect appears when, at a particular value of z the second order contribution begins to 
be comparable to the previous orders and the expansion is no longer valid. Note that in the 
range of low z, which is the most experimentally interesting, the approximation from Table 
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II is thus proved to be acceptable. Similar properties can be observed in the expansion of 
the function R 6 (x, r, (, A). 

For neutral ZA-mesons the situation is very similar to the one with il-mesons. In order to 
study the violation of (12411 . it is necessary to consider the dependence on z of the functions 
R 56 (x, r , C, A). However, unlike the K -meson case, the parameters r and ( for ZA-mesons 
are not well fixed from experiment. The evolution of the set of parameters r and (, which 
violate (1241) . with z (or ct ) is shown in FIG. [3j Gray areas correspond to the function 
Rs(x, r, £, A). Hatched areas correspond to Re(a:, r, (, A). The area of experimentally 
allowed values of the parameters r and ( is contained within the rectangle. At t — 0 the 
areas do not intersect, but have only a common point at r = 1 and ( = 0°. As t —> +oo 
both areas shrink to the point r = 1 and ( = 0°, corresponding to the results of Table II. 

From FIG. [3] one can see that with r > 1 in the limit t —> 0, inequality (|24|1 is only 

violated for the function R 5 (;r, r, (, A). For r < 1 in the limit t —> 0 the violation only 

happens for the function Re (a;, r, £, A). This statement is illustrated in FIG. 0] For the 

function Rs(x, r, (, A), r is set to 1.1. For the function R 6 (x, r, (, A), r is set to 0.9. The 

value of ( in both cases is set to —10°. In analogy with the case of JF-mesons, there is a 
restoration of the inequalities (124|) . but with a higher value of z, ~ 280, which is beyond 
experimental reach. Maximal violation of (124|) also happens at £ ~ 150. From FIG. [5] one 
can see that in the experimentally allowed area z < 3 the violation of (l24ll is less than 10%. 

The difference in behaviour of the functions Rs^^, r, £, A) for K— and ZA-mesons is 
linked to the value of the ratio |Ar|/r, which sets the scale of the horizontal axis. For 



2, while for ZA-mesons this parameter is smaller by almost two 


K mesons, 


K 



As was pointed out above, the study of the violation of (12T11 for ZA s -meson systems requires 
the functions R 7 ; s(^, r, A). In FIG. [6] we show how the areas of violation of (1241) depend 
on z or ct for the functions R. 7 (a;, r, £, A) (gray area) and Rs(x, r, £, A) (hatched area). The 
vertical band shows the experimentally allowed values of r and (. For t = 0 the areas have 
a single point of intersection, r = 1 and ( = 180°. As t —> +oo they shrink to a point at 
(1, 180°). Unlike the case for ZA-mesons, the areas of violation for the ZA s -mesons do not 
evolve monotonical. This is due to the oscillations of ZA s -mesons, which play an important 
role here. 
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FIG. 1: Functions R.s.etx, r, (, A) for neutral A"-mesons (both functions are almost juxtaposed due 
to the high accuracy of the C-P-violation parameter e). The scale at the top corresponds to the 
variable ct (mm); the bottom scale - to the time in units of the average lifetime z = (r^Tr^) t/2 = 
r f, where t is kaon rest frame time. 


ct (mm) 



FIG. 2: Functions r, (, A) for neutral AT-mesons for z <3. 







It is experimentally established that for P s -mesons, r > 1. Hence for z —> 0 only the 
function R 7 (a;, r, (, A) violates the inequality flMll . However in FIG. [8] one observes that 
for z > r/(2A M) ~ 1/2 the function R 8 (x, r, (, A) also begins to violate (I23|) . For the 
numerical simulation, the following values of the parameters were used: r = 1.004 and 
( = 185°. The maximum violation of is reached in the area z ~ 20. At z ~ 40 the 
inequalities are not violated, as in D-meson systems. Due to the high value of z this effect 



0.13, the corresponding values of z are 


intermediate between the ones for K- and D-mesons. 

Conclusions 

Using the oscillations of neutral pseudoscalar mesons we demonstrate the advantages of 
the time-dependent Wigner inequality (J5]) over the static inequality (0J. Eight new time- 
dependent inequalities (124T) were obtained. They can be violated by proper choices of AT 
and i/p for K-, D- and R s -mesons. Relaxation of the obtained inequalities at high values of 
the variable z, is found. This effect is governed explicitly by the CP-violation parameters 
of the considered systems. The inequalities (1241) may be tested at contemporary high-energy 
physics experiments 

V. ACKNOWLEDGEMENTS 

The authors would like to express deep gratitude to Dr. S. Baranov (Lebedev Physical 
Institute, Russia) for numerous discussions related to tests of the Bell inequalities in particle 
physics; to Dr. A. Grinbaum (CEA-Saclay, SPEC/LARSIM, France) for educational chats 
on the foundations of quantum theory; to Prof. S. Seidel (University of New Mexico, USA) 
for help with preparation of the paper. 

We thank for support the Russian Ministry of Education and Science (grant 


N14.610.21.0002, ID RFMEFI61014X0002) and the Russian program of support of the lead¬ 
ing scientific schools (grant SS-3042.2014.2). 


15 



Appendix A: Probabilities required for the time-dependent Wigner inequalities 


In this appendix we summarize all the probabilities that are required to obtain the static 
(J3D and time-dependent © Wigner inequalities for the correlated systems of neutral pseu¬ 
doscalar mesons. 

In the framework of quantum theory, using the normalization condition (J6]) and the 
initial condition ([7]), one can obtain the following expressions for the time-independent 
probabilities: 
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P-?| 2 = ^(1 - cosC sin(2/3)); 
p-g| 2 = i(l - cosC sin(2/3)); 
P + q\ 2 = ^(1 + cosC sin(2/3)); 
P\ 2 = ^ cos 2 /3; 

,| 2 = i S m 2 ft 
p| 2 = i cos 2 ft 
9 | 2 = ism 2 /3. 


To obtain the inequality (J5J) for correlated pairs of mesons MM, we need the values of 
the following time-dependent probabilities (to = 0 below): 

1 l q p 
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~ 4 6 


- 2 rt 
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Appendix B: Functions Rn 

The functions Rn(A, r, £, A) depend on the dimensionless parameters x = ATf, A = 
AM/ Ar , the absolute value r, and the phase C of i/p. In the most general way these functions 
can be written as: 

Rn(z, r, C, A) = fW + foP cos(2Aa;) + / 0 ( f ) sin(2Ax) + / 1 ( f ) sh(x) + / x ( f } ch(x) + 

+ fSc sh (|) cos(Aa;) + ch cos(\x) + ch sin(Ax), 

where we explicitly show the dependence of Rn on x and A. The coefficients /^(r, () are: 
for N = 5 
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for N = 6 
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FIG. 3: Areas of violation of the inequalities (|24l) for D-mesons for the functions Rs^x, r, (, A) 
in the parameter plane of r and ( (( is measured in degrees) depending on z or ct. The gray areas 
correspond to the violation of the function Rs(x, r, (, A), while the hatched areas correspond to 
the function Re(x, r, (, A). The experimentally allowed area of r and £ is contained within the 
rectangle. 
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FIG. 4: Functions Rs^Oe, r, (, A) for neutral D-mesons. The scale at the top corresponds to ct 
(mm); the bottom scale corresponds to the time in units of the average lifetime z = (r#+T l) t/2 = 
r t, where t is calculated in the D-meson rest frame. One can see that with the proper choice of 
the functions Rat for r > 1 and r < 1 the time-dependent Wigner inequalities (|24l) are violated in 
the whole experimentally accessible range of z. 
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FIG. 5: Functions Rs^x, r, (, A) for neutral D-mesons in the area z < 3. The ct here is measured 
in microns (/im). 
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FIG. 6: Areas of the violation of (1241) for B s -mesons for the functions r, (, A) in the r-( 

plane (£ is measured in degrees). The gray areas correspond to the function R 7 (x, r, (, A) the 
hatched areas correspond to the function Rs(x, r, A). The vertical band corresponds to the 
experimentally allowed area of r and (. 
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FIG. 7: Functions Ry^z, r, (, A) for B^-mesons. The scale at the top corresponds to the ct (mm), 
while the bottom scale corresponds to time in units of the average lifetime z = (Th + Tl) t/2 = Tt, 
where t is calculated in the B s -meson rest frame. One can see that the time-dependent inequalities 
(1241) 1 are violated (taking the proper R^r for r > 1 and r < 1) in almost all of the experimentally 
allowed range of 2 . 
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FIG. 8: Functions Ry^x, r, Ci f° r -mesons in the range z < 3, which is the most experimen¬ 
tally interesting. Both functions are almost consistent with one, while the Rs slightly exceeds one 
at z —>• 0. Unlike FIG. [71 here the ct is in microns (/xm). 
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